HD 181068: A Red Giant in a Triply-Eclip 
Compact Hierarchical Triple System 
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Hierarchical triple systems comprise a close binary and a more distant compo- 
nent. They are important for testing theories of star formation and of stellar 
evolution in the presence of nearby companions. We obtained 218 days of Ke- 
pler photometry of HD 181068 (magnitude of 7.1), supplemented by ground- 
based spectroscopy and interferometry, which show it to be a hierarchical 
triple with two types of mutual eclipses. The primary is a red giant that is in a 
45-day orbit with a pair of red dwarfs in a close 0.9-day orbit. The red giant 
shows evidence for tidally-induced oscillations that are driven by the orbital 
motion of the close pair. HD 181068 is an ideal target for studies of dynamical 
evolution and testing tidal friction theories in hierarchical triple systems. 

The Kepler space mission is designed to observe continuously more than 10^ stars, with the 
ultimate goal of detecting a sizeable sample of Earth-like planets around main-sequence stars 
(7). We obtained 218 days of Kepler photometry {2-4) of HD 181068, a star with magnitude 
V = 7.1 and a distance of about 250 pc. It has been previously identified as a single-lined 
spectroscopic binary (5) but there have been no reports of eclipses. 

The data were obtained in long-cadence (LC) mode (one point every 29.4 minutes) over 218 
days using Quarters 1, 2 and 3. Our observations reveal a very distinctive light curve. It shows 
eclipses every ~22.7 days and slow variations in the upper envelope (in Fig. that are likely 
caused by ellipsoidal distortion of the primary component. There are also very regular and much 
narrower eclipses (Fig.flt). These minima have alternating depths, corresponding to a close pair 
(B and C), with an orbital period of ~0.9 d. The 22.7-d eclipses all have similar depths, but there 
are subtle differences between consecutive minima. Radial velocity observations (Supporting 
Online Material, SOM) show that the true orbital period of the BC pair around the A component 
is 45.5 d. The narrow 0.9-d eclipses essentially disappear during both types of the deep minima, 
implying that the three stars have very similar surface brightnesses, so that when the BC pair 
is in front of A, their mutual eclipses do not change the total amount of light coming from the 
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Figure 1: Kepler-hand light curve of HD 181068 from observations in long-cadence mode: (a) 
the full 218 days; (b) a 28-d segment showing two consecutive deep minima, (c-f) close-ups 
of two secondary minima and two primary minima of the 45.5 d eclipses. The dashed and 
dotted lines mark the primary and the secondary minima of the 0.9 d eclipses, respectively. The 
discontinuities in the top panel correspond to the telescope rolls at the end of each quarter. 
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Table 1: Orbital elements for the wider system derived from the A component's radial velocity 
curve (Fig.lS3l). 

Element Value 

Pa-bc 45.5178 d (fixed) 

TMini 2455454.573 ± 0.095 (fixed) 

K2 37.195 ±0.053 km s-i 

v^ 6.993 ±0.011 km 

62 0.0 (fixed) 

/(m) 0.24 ±0.02 Mo 



system (in accordance with the nearly-equal depths of the two deep minima). When the BC 
pair is in front of A, the EC's secondary eclipses appear as tiny brightenings (Fig. [Hi and f), 
showing that the surface brightness of B is almost equal to that of A, while C is a bit fainter, so 
that its disappearance behind B allows the extra light from A to reach us. 

The observed changes in the eclipses of the BC pair and the radial velocity variations of the 
A component confirm that the A and BC systems are physically associated and not a chance 
alignment. Their periods are Pbc=0-90567(2) d and Pa-bc = 45.5178(20) d. Given the shallow 
depths of the eclipses, star A must be the most luminous object in the system. In addition to 
the eclipses, there are brightness fluctuations during the long period minima which imply that 
component A is also an intrinsic variable star with a mean cycle length close to half the shorter 
orbital period, possibly indicating tidally-induced oscillations. 

In addition, there were several flare-like events in the light curve that usually lasted about 
6-8 hours. We checked the Kepler Data Release Notes (7) for documented instrumental effects 
in the vicinity of the 'flares', but found none. Moreover, almost all flares appear right after the 
shallower minimum of the BC pair, suggesting that this activity might be related to the close 
pair. 

We looked for optically resolved companion(s) with a 1-m telescope (SOM), but found 
none. We also obtained 41 high-resolution optical spectra to measure the orbital reflex motion 
of the A component (SOM). The orbital parameters for the wider system (Table 1) reveal that 
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Figure 2: Squared visibility versus spatial frequency from PAVO on CHARA. Grey points show 
all collected measurements, and black symbols the average of each scan over all wavelength 
channels. Each symbol type corresponds to a different night of observations. The solid line is 
the best fitting model. The inset shows a close-up of the observations at the longer baselines. 

star A revolves on a circular orbit, which has an orbital period twice the separation of the two 
consecutive flat-bottomed minima in the light curve (SOM). Long-baseline interferometry using 
the PAVO beam combiner (Precision Astronomical Visible Observations, (S)) at the CHARA 
Array (Center for High Angular Resolution Astronomy, (9)) show that the angular diameter 
of HD 181068 A, corrected for limb darkening (10) is 9ld = 0.461 ± 0.011 milli-arcsecond 
(Fig. [21). 

Combining the measured angular diameter with the Hipparcos parallax of 4.0 ± 0.4 mas 
(11), we find the linear radius of the primary component to be i? = 12.4 ± 1.3i?0. Using 
the spectroscopically determined Teg = 5100 ± 200 K, this a luminosity of L = 93 ± 19Lq. 
This value matches that found from the Hipparcos parallax and the apparent magnitude. We 
also estimated the absolute magnitude of HD 181068 A based on the Wilson-Bappu effect (12), 
which correlates the width of the chromospheric Ca II K emission line at 3934 A with the V- 
band absolute magnitude. Using the latest calibration (S5), the measured width of the emission 
core Wo=12.S km/s implies an absolute brightness of My = —0.3 mag, which matches the 
parallax and the interferometric results. 
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We estimated the mass of HD 181068 A by comparing the effective temperature and lumi- 
nosity with evolutionary tracks from the BASTl database {14). We obtained Ma ~ 3.O±O.4M0, 
corresponding to a red giant, possibly in the He-core burning phase of its evolution {15). The 
full spectral energy distribution, constructed using all published broad-band optical magnitudes 
and infrared flux values, does not show any excess in comparison to a 5200 K photospheric 
model, indicating that no detectable circumstellar dust, possibly from mass-loss processes on 
the red giant branch, is present. 

We have constrained the parameters of the BC pair by modelling the short-period eclipses 
in the Kepler band using the JKTEBOP code {16, 17). The ratio of the radii of the B and C 
components is poorly constrained at present, partly because of the low sampling rate of the 
Kepler long-cadence data. The A component contributes 99.29% of the system light in the 
Kepler passband, and the BC pair contribute 0.44% and 0.27%, respectively. Taking the V- 
band absolute magnitude of HD 181068 A to be My (A) = —0.3 and assuming that our results 
for the Kepler passband are representative of the F-band, we find My (B) = 5.6 and My (C) = 
6.1. Such absolute magnitudes indicate spectral types of G8 V and Kl V for stars B and C, 
respectively {18). Since we do not have independent measurements of Tefr for the BC pair, we 
can only estimate their masses based on their spectral types. This indicates that their masses are 
about 0.7±0.1 Mq each (SOM). 

One puzzling feature of the system is the short-period fluctuations that have the largest 
amplitudes when the BC pair is behind star A, while remaining apparent with a slowly changing 
amplitude in all the other phases of the wide orbit. We have investigated this variability of 
HD 181068 A with a detailed frequency analysis and a comparison to other red giant stars that 
have similar properties (SOM). The frequency content of the light curve suggests an intimate 
link to tidal effects in the triple system, with the first four dominant peaks in the power spectrum 
identifiable as simple linear combinations of the two orbital frequencies. On the other hand, 
solar-like oscillations (meaning those excited by near-surface convection, as in the Sun but also 
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observed in red giants) that are expected to produce an equidistant series of peaks in the power 
spectrum, are not visible, even though all stars with similar parameters in the Kepler database 
do show clear evidence of these oscillations, hi other words, the convectively driven solar-like 
oscillations that we would expect to see in a giant of this type seem to have been suppressed. 

In a recent compilation of 724 triple stars (19), there is only one system with an outer or- 
bital period shorter than that of HD 181068 (A Tau, for which Pout=33.03 d). (20) has reported 
the discovery of KOI-126 with similarly short outer period (Pout=33.92 d). Extremely com- 
pact hierarchical triple systems form a very small minority of hierarchical triplets, with only 
7 of the catalogized 724 systems having outer periods shorter than 150 days. Furthermore, 
HD 181068 and KOI-126 have the highest outer mass ratios (~2.1 and 3, respectively, defined 
as uia/ {mB + mc) among the known systems. In 97% of the known hierarchical triplets be- 
fore the Kepler era, the mass of the close binary exceeded that of the wider companion, and 
even the larger outer mass ratio remained under 1.5. It is not yet clear if this rarity of such 
systems is caused by an observational selection effect or has an underlying stellar evolutionary 
or dynamical explanation. 

Its properties make HD 181068 an ideal target for dynamical evolutionary studies, and for 
testing tidal friction theories. Because of its compactness and its massive primary, we can expect 
short-term orbital element variations on two different time-scales of 46 days (i.e. with period of 
^A-Bc)> and approximately 6 years (^a-bc/^bc)> the time-scale of the classical apsidal motion 
and nodal regression (27), which for the triply eclipsing nature, could be observed relatively 
easily. 
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Sections to the Supporting Online Material 
A Observations and methods 

A.l Lucky Imaging 

We checked HD 181068 for resolved optical companions with lucky imaging on the Im RCC 
telescope of the Konkoly Observatory. For this, we took over 100,000 short-exposure frames on 
2010 June 28/29 and June 29/30, using an Andor IXon^^+888 EMCCD, with exposure times of 
30-61 ms in U BV{RI)q filters. The median seeing was about 1.6". In each filter we obtained 
10,000-30,000 frames, from which the best 0.3% was selected and combined. The resulting 
images show clean single-star profiles with typical FWHM of 0.9" in U, 0.64" in V and 0.45" 
in Iq. For the latter, we have also determined contrast limits at representative separations. An 
analysis of simulated artificial stars with a wide range of brightnesses and separations resulted 
in the following upper limits to the magnitude difference of a hypothetic optical companion: 
0.4" - 2.1 mag; 0.8" - 2.9 mag;1.2" - 5 mag; 2.0" - 6.4 mag. 

A.2 Spectroscopy 

To measure the orbital motion of HD 181068 A, we acquired optical spectra at four different 
observatories. We obtained 41 spectra in total, as follows: 6 spectra with the FIES spectrograph 
at the Nordic Optical Telescope (NOT; resolution 47000, wavelength range 3623-7270 A); 14 
spectra at the Dominion Astrophysical Observatory (DAO; resolution 10000 and wavelength 
range 4300-4556 A); 16 spectra with the 2-m telescope at the Thiiringer Landesstemwarte 
(TLS) in Tautenburg (resolution 66000, wavelength range 4700-7400 A); and 5 spectra at the 
McDonald Observatory (McD) using the 2.7m telescope and the Robert G. TuU coude spectro- 
graph (resolution 60000 and wavelength range 3700-10000 A). 

From the light curve, we know that HD 18 1068 A contributes almost all the light of the 
triple system. We fitted theoretical template spectra from the library of {SI) to a NOT spectrum 
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and determined the following parameters: T^s = 5100 ± 200 K, log^ = 2.8 ± 0.3, [M/H] = 
—0.6 ± 0.3 and v sini = 14km s^^. These values are all in good agreement with those of (52). 
As a check, we can also use Stromgren photometry from (53): V = 7.091, b — y = 0.586, 
mi = 0.296 and Ci = 0.399. From these quantities and the calibration of {S4), we find: 
Teff = 5200 K, E{b - y) = 0.087 and {b - y)o = 0.499. These numbers are consistent with the 
spectroscopic results and confirm that HD 181068 A is a G-type giant star. 

We have used one high-quality NOT spectrum to measure the width of the Ca II K emission 
line at 3934 A. This has been done interactively with the IRAF task SPLOT after converting the 
wavelength scale to Doppler velocities around the core of the emission. The measured value is 
Wo = 72.8 km s^^, which was then used to calculate the V^-band absolute magnitude via the 
calibration of the Wilson-Bappu effect by (55). 

Radial velocities were determined using the IRAF task FXCOR. The template for all but the 
TLS spectra was selected from (SI), with closely matching parameters, which ensured that no 
systematic errors were introduced by spectral template mismatch. The 16 TLS RVs have been 
determined in a first step from cross-correlation with the mean spectrum that was iteratively built 
from the single spectra by shift-and-add according to the measured RVs. This mean spectrum 
was then analysed using the program LLmodels (5(5) to compute a grid of stellar atmosphere 
models and the program Synth V (57) to compute the synthetic spectra. We found Teg = 
5300 ± lOOK, logg = 2.8 ± 0.2 dex, [M/H] = -0.2 ± 0.1, and w sinz = 14 ± 1 km s^^. 
Finally, we used the best-fit synthetic spectrum as a template for cross-correlation to determine 
the RVs of the TLS spectra on an absolute scale. Depending on the instrument and the spectra, 
the velocities are accurate to ±0.5-2 km s^^. The observed radial velocities and their deviations 
from the orbital fit are listed in Table [S2l 

The orbital solution (see Fig. [S3l) was calculated by the method of differential corrections. 
We omitted the DAO RVs because they show a much larger scatter around the calculated orbital 
curve, despite being observed during the same epoch as the other instruments. The eccentricity 
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Table S2: The list of radial velocity measurements. The typical uncertainty is ±1 km s ^. 
Sources: 1=TLS, 2=N0T, 3=MCD, 4=DAQ 
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Figure S3: Measured RVs versus the orbital phase (45.5 d) for the TLS (circles), NOT (trian- 
gles), McDonald (McD, crosses), and DAO (pluses) observations. The vertical bar in the lower 
left comer shows the size of the representative ±1 km s^^ uncertainty. 

from the fit was e = 0.022 ± 0.023, which is consistent with zero, and in the final solution 
we set e = because its inclusion as a free parameter did not improve the solution. We also 
fixed the orbital period to that obtained from the light-curve fitting (a free search gives a slightly 
different value but does not improve the quality of the solution). Table [T] lists the derived orbital 
elements. In the corresponding solution, we corrected the NOT RVs by —0.35 km s^^ and the 
McDonald RVs by —0.70 km s^ with respect to the TLS RVs. This correction minimized the 
rms to 214 m s^\ and its influence on the derived elements was only marginal. 

A.3 Interferometry 

Interferometric observations were performed on three nights in 2010 July, using two different 
baselines (156.3 m and 248.1 m) of the CHARA Array. All the observations were performed 
outside the long-period eclipses, meaning that some flux from the BC pair was present during 
all observations. With an eclipse depth of only 1%, however, the companions are much fainter 
than the primary and are negligible in our analysis. The raw data were reduced using the RWO 
data analysis pipeline. 

We obtained 7 scans of HD181068 over 23 wavelength channels spanning from 0.65 to 0.8 
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mum each, yielding a total of 161 V2 measurements. We observed six stars of spectral type AO 
to calibrate the raw visibilities. All stars were located within 5° in the sky and their estimated 
diameters were at least a factor of 2 smaller than HD 181068A. Inspection of the data revealed 
that three of the calibrators (HD179395, HD182487 and HD181521) are potential binaries and 
could not be used for calibration. The remaining three stars used to calibrate our data were 
HD179733, HD180138 and HD184787. To determine the angular diameter of HD181068 A a 
limb-darkened disc model (S8) was fitted to the calibrated visibilities. The corresponding linear 
limb-darkening coefficient was determined by interpolating the spectroscopically determined 
values of log g and T^s in the grid of (59), yielding /i = 0.63 ± 0.02. The resulted angular 
diameter of HD181068 A is Old = 0.461 ± 0.011 milli-arcsecond. 

The uncertainty on the diameter was estimated using 40000 Monte-Carlo simulations as 
follows: For each simulation, we drew realizations corresponding to Gaussian distributions for 
the calibrator angular diameters, limb-darkening coefficient and wavelength channels, with as- 
sumed standard deviations of 5%, 3% and 0.5%. With these parameters we then calibrated the 
raw visibility measurements and fitted an angular diameter to the calibrated data using least- 
squares minimization. To account for random errors, including correlations between wave- 
length channels, we perturbed this fit by adding random numbers generated from the empirical 
covariance matrix of the data and then repeated the fit (SIO). The final uncertainty was taken 
as the standard deviation of the resulting total distribution, scaled by square root of the reduced 

value (2.5) as determined from the fit to the original data. The measurements of the 7 scans 
averaged over 23 wavelength channels are listed in Table [S3l 

A.4 HD 181068 B and C 

We have constrained the parameters of the BC pair by modelling the short-period eclipses in the 
Kepler band. First, we removed the long-term variations of the uneclipsed brightness in the light 
curve by fitting spline function polynomials and removing data obtained during the long-period 



15 



Table S3: Calibrated interferometric measurements of HD181068 averaged over 23 wavelength 
channels for each scan. The full list of measurements are available from the authors on request. 



JD-2451545 


Projected Baseline (m) 


Spatial Frequency (xlO^rad ^) 


avg(V2) 




3837.973 


144.064 


2.010901 


0.648 


0.040 


3838.000 


140.892 


1.966629 


0.593 


0.043 


3852.920 


247.766 


3.458418 


0.270 


0.050 


3852.930 


247.654 


3.456854 


0.252 


0.053 


3853.929 


247.640 


3.456666 


0.252 


0.044 


3853.944 


247.378 


3.452999 


0.210 


0.057 


3853.969 


246.584 


3.441913 


0.204 


0.054 



eclipses. In this way, the light from star A was assigned to be the 'third light' component. The 
parameters of the BC pair were deduced by modelling the short-period eclipses in the Kepler 
band using the JKTEBOP code {Sll, S12). A preliminary fit was performed, allowing a few 
outlying data points to be identified and removed. A detailed fit was then made, using numerical 
integration to account for the 30-minute long duration of individual observations. Uncertainties 
in the parameters of the fit were calculated using 1000 Monte Carlo simulations (513). In each 
simulation, a synthetic dataset was created by evaluating the best-fitting model at the observed 
times and adding Gaussian noise. This was then fitted in the same way as the real data, starting 
from initial parameter values which were perturbed versions of the best- fitting parameter values. 
The error estimate for each fitted parameter was then evaluated by taking the inner 68.3% of 
all the values of the parameter found from the synthetic datasets. The resulting photometric 
parameters are given in Table IS4l 

We found the ratio of the radii of the B and C components to be poorly constrained at 
present, partly due to the low sampling rate of the Kepler long-cadence data. The A component 
contributes 99.29% of the system light in the Kepler passband, and the BC pair contribute 
0.44% and 0.27%, respectively. Taking the V-band absolute magnitude of HD 181068 A to be 
Mv(A) = —0.3 and assuming that our results for the Kepler passband are representative of the 
y-band, we find Mi/(B) = 5.6 and My (C) = 6.1. Such absolute magnitudes indicate spectral 
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Table S4: Photometric parameters obtained for the short-period eclipses, using standard sym- 
bols in the study of eclipsing binary systems. La,Lb and Lq are all expressed as fractions of 
the total light of the system in the Kepler passband. 



Parameter 


Best fit 


Uncertainty 


Pbc (d) 


0.9056770 


0.0000026 


TMinI (BJD) 


2455051.23625 


0.00020 


ii (degrees) 


87.7 


1.6 


{Rb + Rc)/aBC 


0.3288 


0.0044 


Rc/Rb 


1.01 


0.13 


Rb/ (^bc 


0.164 


0.011 


-Rc/«BC 


0.165 


0.011 


La 


0.9929 


0.0006 


Lb 


0.0044 


0.0007 


Lc 


0.0027 


0.0004 



types of G8 V and Kl V for stars B and C, respectively. 

We can only estimate masses of the BC pair based on the spectral type because there is 
no independent Teg measurement of them. Fig. [S4] shows a Hertzsprung-Russel Diagram with 
BASTI evolutionary tracks (S14) of different masses for [M/H] = —0.25 (as suggested from 
the spectroscopy), as well as the most likely locations of the components of the triple system. 
The position of the BC components was derived based on the spectral types stated above, as- 
suming an error of 200 K in T^s and 0.2 mag in absolute magnitude, respectively. 

B The variability of HD 181068 A 

The light curve (see Fig. [1]) shows slow variations with the same timescale as the long-period 
eclipses, which presumably arise from ellipsoidal distortion of the primary. We also see faster 
oscillations with the same timescale as the orbital period of the BC pair. These are visible both 
outside and during the eclipses and are less obvious to interpret. 

To investigate this further. Fig. ISSb shows the amplitude spectrum of the light curve after 
first removing observations made during both the long- and short-period eclipses. This proce- 
dure left a light curve with a duty cycle slightly above 60%, and it introduced alias peaks in the 

17 



3.9 3.8 3.7 3.6 

log T,„/K 

Figure S4: A Hertzsprung-Russell Diagram using the BASTI tracks {S14) for \M / H] = —0.25. 
It shows the locations of components A, B and C, marked by the filled square, diamond and 
circle, respectively. 

amplitude spectrum at the multiples of the orbital frequencies. The strongest peak in the spec- 
trum occurs at 25 /xHz, corresponding to half the orbital period of the BC binary. As mentioned, 
this periodicity is clearly visible in the light curve. 

To search for other frequencies we used iterative sine-wave fitting (prewhitening) with Pe- 
riod04 (Si 5). In five steps, we measured and identified the following frequencies: 



• 


/l 


= 24.54 /iHz 


2(/gijQ[.t 2/iQng 


• 


/2 


= 25.05 /iHz 


2(/ghort /long) 


• 


/3 


= 25.56 /iHz 


2/gjjQj.t 


• 


h 


= 51.12 /iHz 




• 


h 


= 12.83 /iHz 


= /short + 1/^obs 



Here, /short = 1/Pbc, /long = I/^a-bc, and Tobs is the time span of observations. Fig.|S5^ 
shows the amplitude spectrum after the slow variations and these five strongest peaks have been 
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subtracted from the time series. 

Could the observed signal in HD 181068 arise from solar-like oscillations? According to the 
parameters derived from the interferometry and spectroscopy, HD 181068 A is a red giant star 
located in the bottom of the red giant branch. Studies of similar red giants with Kepler {S16, SI 7) 
show that essentially all stars in this region of the H-R diagram exhibit solar-like oscillations 
that appear as a broad power excess centred at a characteristic frequency z/max and composed of 
a regularly spaced series of peaks. Using the scaling relation of (SI 8) with the derived values 
of mass, radius and effective temperature, we estimate z/max = 64 ± 16 /iHz for HD 181068 A, 
while the expected amplitude is about 80 ppm. 

For comparison, Fig.lSSb shows the amplitude spectrum of a typical red giant with z/max ~ 
80 /xHz (KIC 12507577). To make the comparison exact, we calculated this amplitude spectrum 
using exactly the same portions of the light curve that were used for HD 181068 in Fig. \S5[ i. 
In KIC 12507577 we see the regular peaks that characterize solar-like oscillations, whereas in 
HD 181068 A we see just a few peaks whose removal (Fig. \S5h ) leaves only a slowly rising 
power distribution. The observed signal in HD 181068 A is clearly not compatible with solar- 
like oscillations. Indeed, the solar-like oscillations that we would expect to see in a giant of this 
type seem to have been suppressed. 

The frequency content of the light curve suggests an intimate link to the orbital frequencies 
in the triple system. We are led to suggest that we are seeing tidally-induced oscillations that 
are driven by the orbital motion of the BC pair. Tidally-induced oscillations have previously 
been reported in a few binary systems {S19, S20, S21), but here the situation is different be- 
cause the period of the oscillations does not correspond to the orbit of the A component, but 
rather to that of the BC pair. While a fuller discussion of this possibility is postponed to a 
future publication, we note that the amplitude of tidally driven oscillations can be simply es- 
timated by assuming that the brightness changes are proportional to the tidal heights given by 
Ra ( ( Mb + Mc ) / Ma ) ( Ra /a)^, where a is the semimaj or axis of the outer binary. Every num- 
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ber put together yield an estimated amplitude of 3 ppt, which is a factor of 10 larger than what 
we observe and a factor of 40 larger than the amplitude of solar-like oscillations. 

C Supplemented data files 

In addition to this document describing the supporting extra material, we also make all the data 
mentioned throughout the paper and SOM available to the general community at the webpage 
of the journal. The full list of datafiles that are stored in a single compressed tarball file is the 
following: 

12 017 62som_sed.txt 
12017 62som_luck.y_image . fits 
12 017 62som_fies (1 . . 6) . fits 
12 017 62som_dao (1 . . 11) . fits 
12 017 62som_tls (1 . . 16) . fits 

1201762som_mcd (1 . .5) .fits 
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Figure S5: The amplitude spectrum of HD 181068 A and the change after five prewhitening 
steps (panels a and b). The inset shows the spectral window function with the same axis scales. 
Panel c: the spectrum of a typical red giant star (KIC 12507577), for which v^ax ~ 80 /xHz. 
The vertical scale in panels b and c is increased by a factor of two. 
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